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ABSTRACT 

We investigate the relevance of the gravitational-wave driven r-mode instability for 
strange stars. We find that the unstable r-modes affect strange stars in a way that 
is quite distinct from the neutron star case. For accreting strange stars we show that 
the onset of r-mode instability does not lead to the thermo-gravitational runaway that 
is likely to occur in neutron stars. Instead, the strange star evolves towards a quasi- 
equilibrium state on a timescale of about a year. This mechanism could thus explain 
the clustering of spin-frequencies inferred from kHz QPO data in Low-mass X-ray 
binaries. For young strange stars we show that the r-mode driven spin-evolution is 
also distinct from the neutron star case. In a young strange star the r-mode undergoes 
short cycles of instability during the first few months. This is followed by a quasi- 
adiabatic phase where the r-mode remains at a small, roughly constant, amplitude for 
thousands of years. Another distinguishing feature from the neutron star case is that 
the r-modes in a strange star never grow to large amplitudes. Our results suggest that 
the r-modes in strange star emit a persistent gravitational-wave signal that should be 
detectable with large-scale interferometers given an observation time of a few months. 
If detected, these signals would provide unique evidence for the existence of strange 
stars, which would put useful constraints on the parameters of QCD. 
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1 INTRODUCTION 



Quite plausibly, the most stable form of matter at supranu- 



clear densities is a conglomerate of deconfined strange, up 
and down quarks. This could have important consequences 
for compact stars, with central densities several time the nu- 
clear satura tion density, that are born following supernova 
explosions (Witten 1984). It has, in fact, been suggested 
that if the strange matter hypothesis is correct then all the 



observed neutron stars ma y be strange stars (Alcock et a] 



198e; Colpi fc Miller 1992). If it could be established that 



strange stars exist in the Universe it would have obvious 
implications for our understanding of pulsars. In addition it 
would put severe constraints on th e QCD parameters [eg. 
the bag constant in the MIT model jFarhi fc Jaffe 1984|)]- It 



is thus important to establish whether observations can dis- 
tinguish a strange star from a neutron star. Unfortunately, 
this turns out to be a delicate problem. 

Strange stars are held together by both the strong nu- 
clear interaction and gravity, and the corresponding equa- 
tion of state is quite acc urately described by uniform den- 
sity models (see Madsen ( 1998b ) for an exhaustive review). 
A consequence of this is that, in contrast to the neutron 
star case, very small strange "dwarfs" (with mass M ^ R^) 



can form. However, for the canonical mass of IAMq grav- 
ity dominates the strong interaction which lea ds to strange 



stars and neutron stars being similar in size (Alcock et a] 



1986 ). In other words, it is not clear that one would be able 
to distinguish between the two cases given observed masses 
and radii. 

One consequence of the extreme stiffness of the strange 
matter equation of state is that strange stars can reach 
higher rotation frequencies than typical neutron stars. How- 
ever, to date the fastest observed pulsar is the 1.56 ms 
PSR1937-(-21 which spins at a rate significantly below the 
breakup limit for most proposed supranuclear equation of 
state. In order to indicate the presence of a spinning strange 
star the period would have to be significantly shorter than 
1 ms. Thus the claimed observation of a 0.5 ms pulsar in the 
remnant of SN1987A caused considerable excitement but, of 
course, it turned out to be flawed. Despite severa l searches 
[see for example (Edwards, Straten & Bailes 2001)] no sub- 
millisecond pulsars have yet been discovered. 

It has been argued that strange stars might be distin- 
guished by the fact that they cool faster than a neutron star 
would, but it is now accepted that the direct Urea mecha- 
nism may be op erating [either by th e proton fraction rising 
above ~ 10 % ( Lattimer et al 1994 ) or by hyperons being 
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present (Page et al 2000)] in the core of a neutron star. This 
would cool the neutron star very rapidly, which would make 
it diffic ult to identify a strange star from cooling data. 



inal r-mode spin-evolution models led to results that accord 
well with the infe rred initial spin period of about 2 ms for 
the Crab pulsar ([Lindblom^^^Owen fc Morsink 199^ ; lAnder 



An other possible indication of a strange star would be sson, Kokkotas fc Schutz 1999| ). However, the discovery of 



thermonuclear flashes from the surface as accreted material 
is converted into strange matter. This mechanism would, 
however, only operate if the strange star was bare. It is gen- 
erally anticipated that strange stars w ill be covered by a 
thin crust of "normal" nuclear matter (Alcock et al 1986). 
This crust is suspended above the strange matter core by 
the ele ctromagnetic field. This moans that the crust densitv 



cannot exceed that of neutron drip (Gleiidonning fc Webei 



1992). As neutrons begin to drip out of the crust nuclei they 
will migrate into the core and be converted into strange 
matter. The presence of such a crust would allow accretion 
to proc eed without necessarily leading to flashes from the 
surface 



The sole observational indicatio n against st range stars 
is provided by the glitching pulsars (Alpar 1987). The stan- 
dard model for explaining the large glitches in for example 
the Vela pulsar relies on the transfer of angular momentum 
between the neutron star crust and a superfluid component. 
It is not easy to see how to construct an analogous model 
for glitches in a strange star. 

Following the discovery that the r-modes in a rotating 
perfect fluid star would be d riven unstable via the emission 



of gravitationa l waves [see ( Andersson fc Kokkotas 2001 
Lindblom 20oTl ) for detailed reviews], Madsen ( |1998at |200C| ) 



suggested that this mechanism would provide the means 
for dist inguishing between neutron stars and strange stars. 



The gis t of Madsen's argument is that the bulk viscosity of 



strange matter is several orders of magnitude stronger than 
that of normal neutron star matter, which essentially shifts 
the window in which the r-modes are unstable to signifi- 
cantly lower temperatures. As a consequence one would not 
expect a newly born strange star to be immediately affected 
by the r-mode instability, while a neutron star might spin 
down significantly in the first few months. An observation 
of a newly born pulsar spinning near the Kepler limit might 
thus provide evidence for the presence of a strange star. 

The realisation that the r-modes in a spinning compact 
star are generically unstable has led to a considerable effort 
aimed at understanding the possible astrophysical relevance 
of this mechanism. In the last two years several crucial issues 



between oscillations in the core fluid and the crust ( Bildsten 


& Ushc 
& Ushc 


mirsky 200C ; 


Andersson et al 2000; 


Lindblom, Owen 


mirsky 2000 


, the role of the magnetic fleld (Spruit 


199S|; |l 


ezzoUa, Lamb & Shapiro 200C; MendeU 200d|), su- 


perfluidity (Lindblom & MendeU 2000l 


Andersson fc Comer 



to 



exist i n the de ep neutron star core (Jones 2001; 



Lindblom 



& Owe [i 2001 ). Of similar importance has been attempts 



to understand the nonlinear saturation of an unstable mode 



via hydrodynamical simulations ( B tergioulas fc Font 2000 
Lindblom, Tohline fc Vallisneri 200o| . For a i nore complete 



set of references, see recent re view articles (Andersson & 
Kokkot [a.s 200l| ; [Lindblom 200l| ) . 

In parallell with the various attempts to implement 
more detailed physics into the r-mode model, there have 
been discussions regarding possible ways that the instability 
may manifest itself in current observational data. The orig- 



the 16 ms PSR 0537-6910 that is thought to have been born 
spinning signiflcantly faster (at 6-9 ms assuming a reason- 
able braking index) indicates that the naive r-mode model 
needs refinement. Furthermore, one can readily show that 
external agents like the fallback of supernova debris and an 
acting magnetic propeller torque are of utmost impo rtance 
in determining the spin of a newly born neutron star (Watts 



Andersson 2001) 



Ado pting a scenario first sugg ested by Papaloizou & 
Pringle ( 1978 ) and Wagoner ( 1984 ), it has been suggested 
that the r-mode insta bilitv may be active in mature ac- 
cretin g neutron stars (Andersson, Kokkotas fc Stergioulas 



199!;). This could provide an explanation for the apparent 
clustering of spin rates (in the range 260-590 Hz) inferred 
from kHz QPO data for Low-Mass X-ray Binaries (LMXB). 
The original ideas were based on the notion that the angu- 
lar momentum lost through gravitational waves radiated by 
the r-modes would balance the accretion torque and prevent 
the star from spinning up. This would lead to a persistent 
gravitat ional- wave signal. Howe ver, as was pointed out by 
Spruit (1999) and Levin (199!;), the onset of instability is 
likely to trigger a thermo-gravitational runaway. This leads 
to a brief period (a few months) of r-mode spin-down fol- 
lowed by a long period (millions of years) of accretion driven 
spin-up. Even though the r-mode instability wo uld still be 
able to cau se the observed clustering of LMXBs ( Andersson 
et al 2000| ), the associated gravitational waves would likely 



not be detected (as the event rate would be too low). 

As we will show in the following, the r-mode instability 
acts rather differently in a strange star. We will argue that 
unstable r-modes in an accreting strange star in a LMXB 
may, in fact, provide a persistent source of gravitational ra- 
diation. We will also consider the evolution of young strange 
stars. We will show how the r-mode driven spin-evolution of 
a strange star differs significantly from the neutron star case. 
Most significantly, we find that the r-modes never grow to 
large amplitudes in a strange star. The veracity of our re- 
sults should be testable with large scale laser interferometers 
such as LIGO, VIRGO and GEO600 that are about to come 
on-line. We thus propose that gravitational- wave astronomy 
might soon be able to establish whether the accreting com- 
pact stars in LMXBs are, indeed, strange and perhaps even 
constrain the parameters of current theoretical models. 

We need to model the evolution of the spin rate, tem- 
perature and r-mode amplitude in an accreting strange star. 
The essential elements of our solution to this problem are 
discussed in sections 2-3, and we present our results for re- 
cycled strange stars in Section 4. In section 5 we discuss 
the evolution of young strange stars. Our conclusions are in 
Section 6. In an Appendix, we discuss possible refinements 
to the model. In particular, we consider the deconfinement 
of the strange star crust as the spin of the star changes and 
how nonlinear contributions to a strong viscous dissipation 
may serve to saturate an unstable r-mode. 
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Table 1. 


Parameterised variables used in the paper. 


Ml.4 


M/1AMq 


mass 


Rio 


R/10 km 


radius 


P-3 


P/1 ms 


rotation period 


Ts 


T/IO* K 


core temperature 


M-8 


Af/lO-SMg/yr 


accretion rate 


620 


e/20 MeV 


energy released per accreted nucleon 


ao.i 


as/0.1 


strong coupling constant 


mioo 


ms/100 MeV 


strange quark mass 


M400 


/id/400 MeV 


down quark chemical potential 


dio 


d/10 kpc 


distance to gravitational wave source 



2 MODELLING UNSTABLE R-MODES IN 
STRANGE STARS 

The aim of this section is to collect results from the litera- 
ture that will enable us to assess the relevance of the r-mode 
instability for, possibly accreting, strange stars. We have at- 
tempted to parameterise all relevant relations in order to 
make the dependence on the stellar (and QCD) parameters 
clear. The various parameterisations are summarised in Ta- 
ble |l[ 



2.1 R-mode estimates 

In the rotating frame of reference the r-modes have fre- 
quency 

2mD. 



(1) 



where Q is the rotation frequency of the star. 

For simplicity we will assume that the star is described 
as an n = 1 polytrope. The main reason for this is that 
the corresponding r-mode timescales have already been sum- 
marised by Andersson & Kokkotas (2001). A more realistic 
strange star equation of state will be stiffer than our model, 
but we do not expect this to affect our results significantly. 

Due to the emission of gravitational waves the I — m = 
2 r-mode (which is expected to lead to the strongest insta- 
bility) grows on a timescale 



to 



-47Mr,lRToP- 



(2) 



(the minus sign indicates that the mode is unstable). We 
include only the contribution from the leading current mul- 
tipole here. 

Viscous dissipation tends to suppress the growth of an 
unstable mode, and we find 



ts 



7.4 X lO'al/fM-l^'Rl'o'-'T^'-' s 



11/3^5/3 



(3) 



for the shear viscosity. For the bulk viscosity (which is due 
to the change in concentration of down and strange quarks 
caused by the mode oscillation) the viscosity coefficient takes 
the form 



(4) 



where the coefficients a and /? are given by Madsen (1992). 
From this we see that the bulk viscosity gets weaker at both 
very low and very high temperatures, cf. Fig. 1 of Mad- 
sen ( 2000 ) . For low temperatures we find that 




Figure 1. The r-mode instability window for a strange star with 
canonical parameters M = IAMq and R = 10 km (solid line). 
For comparison we show the corresponding instability results for 
normal npe fluid and neutron stars with a crust (dashed curves). 
We compare the instability window to i) the inferred spin-periods 
for accreting stars in Low-mass X-ray binaries [shaded box], and 
ii) the fastest known millisecond pulsars (for which observational 
upper limits on the temperature are available) [horisontal lines]. 



Meanwhile, at high temperatures we cannot readily write 
down a parameterised expression for the bulk viscosity 
timescale. The relevant timescale has to be calculated nu- 
merically for each given stellar model. However, this has lit- 
tle impact on our present study since the star would likely 
not remain in the high temperature part of the instability 
window long enough for an unstable mode to grow signifi- 
cantly. Hence, we do not consider the high-temperature case 
(above say 10^ K) further in this paper. 

Given the above estimates one can readily construct a 
critical curve in the fl — T plane above which the r-mode 
will be unstable. The mode is unstable whenever 



t, 



+ '-f 



.low 



-2 -4 

rnjoo s 



(5) 



(6) 

where ig,6,s are the timescales given in Eqns. (^, ^ and 
(^), is negative. If we assume canonical values for the various 
parameters we obtain the results shown in Figure ^ [cf. Fig. 1 
of Madsen ( |2000[ )]. 

From this figure we immediately see that the r-mode 
instability would not become active (apart from a brief ini- 
tial period expected to last a few tens of seconds when the 
star is significantly hotter than 10^ K) until the star has 
cooled to a temperature of a few times 10* K. A young 
strange star would therefore not be spun down significantly 
by the r-mode instability during its first few months of ex- 
istence. But after a year or so the star will enter the "low- 
temperature" instability window and undergo a phase of r- 
mode spin-down. 



2.2 Temperature evolution 

Given that the viscous damping timescales are strongly de- 
pendent on the temperature we obviously need to model 
the thermal evolution of the star. To do this we first of all 
assume that the stellar core is isothermal. 

For a strange star, one can e stimate that the (inte- 
grated) heat capacity is (Usov 2001) 



4 N. Andersson, D.I. Jones and K.D. Kokkotas 



Cv ~ 1.5 X 10=^«Mi'//i?ior8 erg/K (7) 

From this we readily infer that the thermal energy of core is 

i^thcrmai = 7.7 X 10*"i?ioMi'//r8" erg. (8) 

At the temperatures that we will consider, the s tar cools 
mainly due to neutrino emission. Pizzochero (1991) gives a 
formula for the corresponding luminosity per unit volume, 
which, when multiplied by the volume of a spherical star, 
gives: 

^neutrino = 3.77 X 10^'' RloTi ergs-\ (9) 
We can define a cooling rate for an isolated star as 



S-Ethc 



(10) 



(with the factor of 2 included since the thermal energy is 
quadratic in the temperature). Combining equations (^) and 
(H) we obtain 



f « 4.1 X 10*Mj^//i?7cf Tg"* s 



(11) 



(where we have assumed that the initial temperature is high 
enough to make its contribution to the formula insignifi- 
cant). From this formula we see that it would take the star 
roughly 10 years to cool down to a temperature of 10** K. 
A canonical strange star spinning at the Kepler limit en- 
ters the main r-mode instability window in Figure |l| at 
r « 6 X 10* K. Given the above cooling timescale, this 
would happen roughly three days after the star was born. 
The star then evolves through the instability phase for the 
next 10^ years or so, until the core temperature has fallen 
below 10^ K. 

If the star is accreting we must also account for heating 
due to the accreted material. In the case of strange stars, the 
main heating source is likely to be due to the convers ion o f 
accreted nucleons into strange matter. Alcock et al. (198f ) 



estimate a release of e = 20 MeV of heat per nucleon due to 
this conversion. This then gives 



accretion 



ne, 



(12) 



where h is the rate as which nucleons enter the strange mat- 
ter core. As we will discuss in the Appendix this is not nec- 
essarily the same as the accretion rate of nucleons onto the 
star. In the case where all accreted nucleons are immediately 
converted we get 



E, 



accretion 



= 1.19 X 10^^M-8e2oergs \ 



(13) 



Finally, a detailed study should also account for the fact 
that the shear viscosity will contribute to the heating of the 
fluid. We estimate the associated heating rate as 

i? 1 1 ^, in45 2 , ,14/9 J3-5/3 r)-2 -5/3rTi-5/3 

^viscosity = 1.1 X 10 aMj^J R-^q' P_^ctQ{ Tg (14) 

where a represents the amplitude of the r-mode (see below) . 

It is easier to evolve the thermal energy -Ethcrmai rather 
than the temperature. The relevant differential equation is: 



Ea 



on 

+ E 



£'neutrino ; 



(15) 



In the absence of an r-mode, the temperature is deter- 
mined by a balance between accretion heating and neutrino 
cooling. Equations (|l^) and (^ then lead to 

T = 8.25 X 10^ [Af-8e2o-Rro^] K. (16) 



For temperatures of this order, shear viscosity dissipation is 
much smaller than bulk viscosity dissipation, so that the r- 
mode instability sets in when \tg\ — ti,. Equations (^ and (^ 
combine to give the corresponding critical rotation period: 



r,rrn-l/2 , ,3/4 -1 



(17) 



It is illuminating to compare these predictions to the 
inferred clustering in the LMXB data. We predict that stars 
with accretion rates in the range 10~^ < M-s < 1 should be 
confined to spin frequencies in the range 300 — 450 Hz. Given 
the many uncertainties in our simple model, this compares 
quite favourably with the 260 — 590 Hz range suggested by 
observations. Conversely, it is clear that the current obser- 
vational data cannot be used to rule out the possibility that 
the LMXBs harbour strange stars. 



3 THE SPIN-EVOLUTION OF ACCRETING 
STRANGE STARS 

In order to investigate the way in which a strange star 
evolves under infiuence of both accretion and an unstable 
r-mode we use a phenomeno logica l spin evolution model sim- 
ilar to that of Owen et al (1998). We assume that we can 



model the star's angular momentum as the sum of the bulk 
angular momentum J and the canonical momentum of the r- 
mode Jc. The total torque on the star J can then be written 



(18) 



where the dots indicate time-derivatives. In the following 
we will write the moment of inertia as J = IME? with 
I = 0.261 for an n = 1 polytrope. 

Note that, since we are modelling the star as an n = 1 
polytrope, accretion leads to a change in the mass M while 
the radius R remains constant. This is a peculiarity associ- 
ated with our particular stellar model, but it has no effect 
on the outcome of our evolutions. We have verified that this 
is the case by considering other polytropes, and also includ- 
ing the effect of the centrifugal ffattening in our model. This 
more general study leads to results that are virtually iden- 
tical to the ones obtained from the model described here so 
we will discuss only this simple case. 

For the I — m — 2 r-mode one can show that the canon- 
ical angular momentum is 



Jc 



2,Q.a^JMR^ 



(19) 



where a is the mode amplitude (as defined by Owen et al) 
and J = 1.635 x 10"^ for an n = 1 polytrope. 

The r-mode is driven by gravitational radiation and 
damped by viscosity. We thus assume that the canonical 
angular momentum evolves according to 



A 

2Jc 



(20) 



(the factor of 2 is included since Jc is proportional to the 
square of the perturbation). 

The total torque on the star is the sum of torques due 
to gravitational radiation and accretion. Throughout the 
main part of the paper we neglect magnetic dipole radia- 
tion. There are two simple reasons for this. Firstly, such 
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radiation is likely to be suppressed in an accretion environ- 
ment, especially since the magnetic field of the LMXBs are 
unlikely to be stronger than ~ 10^ G. Secondly, in the case 
of young isolated strange stars one can argue that the main 
r-mode spindown torque will dominate the magnetic dipole 
torque. We will comment on this further in Section 7. It is, of 
course, possible that the interplay between the r-mode and 
the interior magnetic field will affect the develop ment o f the 
instability (as sugge sted by, for example Spruit (1999) and 



RezzoUa et al ( 2000 ) ) . We do not account for this possibility 
in our study. 

The torque from the gravitational wave emission from 
the I — m — 2 current multipole is 



Ja 



while we assume accretion to lead to a torque 



Ja = MVGMR 



(21) 



(22) 



This expression is, of course, likely to be a serious simplifica- 
tion of the true accretion torque, but it should be sufficient 
for the rather qualitative considerations of the present pa- 
per. 

Equations (p^-(|2^) combine to give equations for the 
evolution of a and SI. 



-+ 1 



3a^7 
2/ 



M / G 

T \MRy 



MQ 
If 



- 6\la — I 1 

I V tsv ^bv 



G 



1/2' 



(23) 



(24) 



4 LMXBS AS A SOURCE FOR PERSISTENT 
GRAVITATIONAL WAVES 

As long as the r-mode remains stable an accreting star will 
spin up, and its core temperature should be well approxi- 
mated by (^. As is clear from Figure hi this means that the 
star will enter the r-mode instability window on the branch 
where bulk viscosity provides the main damping mechanism. 
The subsequent evolution will therefore be manifestly differ- 
ent from the neutron star case where the star would undergo 
a thermo-gravitational runaway once it reached the instabil- 
ity limit. We shall now show that in the case of strange stars 
the subsequent evolution may well proceed extremely slowly 
(on the accretion timescale) and lead to emission of persis- 
tent gravitational waves. 

As soon as the star enters the instability window the r- 
mode will grow and lead to both shear viscosity heating and 
a spin-down torque. Intuitively, one would expect the star 
to continue spinning up until the r-mode amplitude becomes 
sufficiently large to balance the accretion torque. However, 
the situation will not be stable unless this balance occurs 
for a star residing on the instability curve, i.e. for which we 
have Q = 0. This situation is well approximated hy \tg\ tb- 
Use equation p^ ) and assume that i) the mode amplitude 
is small (such that 1 >> io? J /2I) and ii) the temperature 
is high enough that the shear viscosity dissipation can be 
neglected. 

In order for the situation to correspond to a quasiequi- 
librium the mode-amplitude required to balance the accre- 
tion torque is equal to: 



The differential equations ( p3| ) and (^J) can be inte- 
grated to give trajectories in the SI — T plane along which 
an accreting strange star would be expected to evolve. 

Note that ( ^ ) indicates that the r-mode is unstable, in 
the sense that a > 0, whenever 

1 M ^ 

7+2n + 2M<° (25) 

This is notably different from (^). The difference is, however, 
easily explained. Eqn ( |25| ) follows since Jc is conserved for 
an isolated star, cf. Ho & Lai ( |2QQ(]| ). Eqn (||) is the appro- 
priate criterion for onset of instability in a star that is spun 
up by accretion. In our evolutions the r-mode instability is 
consequently active whenever ( |25| ) holds. 

Equations (|2^) and (ji^) only apply as long as a re- 
mains "small". As an unstable r-mode grows exponentially 
one would expect the equations to remain valid only for a 
very limited period of time. Intuitively, one would expect 
to growth of the mode to be halted once nonlinear effects 
become relevant. Recent attempts to model the evolution of 
large amplitude unstable modes in numerical hydrodynam- 
ical simulations suggest that the r-modes ma y not saturate 
until a ha s grown to values of order unity (Btergioulas & 
Font 200; [Lindblom, Tohline fc Vallisneri 200q ). As we will 
demonstrate below, the r-modes are unlikely to ever grow 
to such amplitudes in a strange star. Thus we assume that 
the r-mode does not saturate due to nonlinear effects, and 
consider equations (^) and (^) to be relevant throughout 
our evolutions. 



; 6.7 X 10~''M- 



11/2 



-7 rp- 

^100 -'8 



■7/2 



(26) 



while the mode amplitude required for the shear viscosity 
heating to prevent the star from cooling down is equal to: 



3.6 X IQ-^U'l'^^^al'^R^o 
[Z.nRloTi - 1.19M-8e2o] 



5/3^-2^2/3 
''''1 nn-t s 



(27) 



In obtaining these relations we have made use of ( |l7| ) to 
express P in terms of T. 

Now we could in principle equate the two expressions 
( pij ) and (^^ and solve for the "equilibrium" temperature. 
However, because this leads to a nonlinear equation it is non- 
trivial. In order to get an algebraic answer we instead assume 
that the star remains close to the equilibrium temperature 
(p^. Linearising the equation with T = Tcq -I- ST it is easy 
to estimate that the required solution corresponds to 



5T_ 

Tea 



0.24 



1 + 942T"'""R 



,25/6 7,43/6 5/3 



10 



0.1 



5 11^-137/36 

^!ooe2oMi,4 



(28) 



From this we can deduce that once the star enters the in- 
stability regime it never evolves far away from the equilib- 
rium temperature given by (p^). In other words, gravita- 
tional waves halt the accretion spin-up once the star reaches 
the r-mode instability curve. 

To complete this picture it is necessary to demonstrate 
that the equilibrium described above is a stable one, i.e. that 
for small departures from the equilibrium values (a, T, P) 
the star returns to the original state. This is straightfor- 
wardly achieved by assembling our three differential equa- 
tions for these quantities, obtaining linearised equations for 
the departure from equilibrium, and writing the time de- 
pendence of this small departure as e™' . The solution to the 
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Figure 2. The evolution of tlie r-mode amplitude (upper frame), 
the temperature (middle frame) and the spin period (bottom 
frame) corresponding to the first 10^" s following the initial onset 
of r-mode instability in an accreting strange star. 




log,oT/lK 

Figure 3. The evolution of a young strange star (with canonical 
parameters M = IAMq and R = 10 km) under the influence of 
the r-mode instability. The thin solid line represents the relevant 
part of the instability window from Fig. |l| The first year of the 
stars evolution following the first onset of instability are shown 
as a thick solid line. 



eigenvalue problem so obtained shows that the configura- 
tion is indeed stable (the real part of w is negative), with 
departures from equilibrium being damped on a timescale 
intermediate between the temperature evolution timescale 
and the gravitational radiation reaction one. This result will 
be of interest should such wandering behaviour ever be seen 
in an accreting system. 

Numerical results supporting these analytic estimates 
are shown in Figure |^. As is readily apparent, after a brief 
interval in which a and T vary sharply on a timescale close to 
the thermal one, the system settles down to a steady state. 
All subsequent evolution occurs on the accretion timescale. 

Having established that the spin-up of an accreting 
strange star will be halted by the r-mode instability we 
can readily assess the detectability of the resultant grav- 
itational waves. Following Ander sson, Kokkotas and Ster- 
gioulas (1999) (see also Bildsten (1998)) we use 



4G 



(-)' 

\LurJ 



(29) 



and E = —uij^cc/m, where uj is the mode frequency mea- 
sured in the inertial frame to give: 



/i « 2.3 X 10~^'^P%^Ml^^Ryo'^M^Jgdio 



(30) 



Since the detectability of the source increases roughly as the 
square-root of the number of observed cycles we find that 
the effective amplitude is 



2tt 



-h ~ 10" 



(31) 



after a few months of observation. In other words: this signal 
should easily be detected by a large scale interferometric 
detector. 



5 THE SPIN-DOWN OF YOUNG STRANGE 
STARS 

We next want to investigate whether the spin-evolution of 
a young strange star is also distinct from that of a nascent 
neutron star. To do this we simply assume that M = in 
the various evolution equations from Section 3 and initialize 
the integration with P = Pk, i.e. the star is spinning at the 
break-up limit when it cools to the temperature at which 
the r-mode first goes unstable. A typical result of this cal- 
culation, corresponding to the first year of the evolution, is 
shown in Figure ^ 

At first sight, the result shown in Figure ^ may be a bit 
surprising. A strange star clearly evolves in a way that is 
rather different from the familiar results for hot young neu- 



tron stars, cf. Owen et al (199S). In particular, the strange 



star never evolves far into the instability window for the 
I = m = 2 r-mode. This is essentially because the main r- 
mode instability phase does not begin until a strange star 
has cooled to a temperature below 10^ K. At that temper- 
ature the cooling time is considerably longer than the spin- 
down time due to an r-mode of sizeable amplitude. This 
means that, once the r-mode begins to grow it efficiently 
spins the star down through the bulk viscosity part of the 
instability curve. Once stable, the mode dies away expo- 
nentially and cooling will again govern the evolution of the 
star. Eventually the mode is again unstable and the cycle 
is repeated. After a few years the system reaches a quasi- 
equilibrium where the mode amplitude evolves on the cool- 
ing time scale, and the star spins down along the bulk- 
viscosity branch on the instability window. That this woul d 
be the case has previously been suggested by Madsen ( 2000 ) . 



This result has considerable repercussions for the r- 
modes as a gravitational-wave source. In Figure ^ we show 
the evolution of the r-mode amplitude corresponding to the 
case shown in figure ^. Here it is interesting to note that 
the mode amplitude varies greatly during the initial phase 
of r-mode spindown. Furthermore, it is clear that the evolu- 
tion never leads to a reaching values of order unity. This is 
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Figure 4. The upper frame shows the evolution of the r-mode 
ampUtude during the first year of evolution of a young strange 
star. It can be seen that the mode never reaches amplitudes of or- 
der unity, and also that the oscillating evolution in the amplitude 
is replaced by a smoother behaviour with small values (~ 10~^) 
of a after a year or so. The middle frame shows the corresponding 
evolution of the gravitational-wave frequency, while the bottom 
frame illustrates the gravitational-wave amplitude for a source in 
the Virgo cluster (at a distance of 15 Mpc). 

intriguing given the results of recent hydrodynamical simu- 
lations that seem to suggest that an unstable r-mode will not 
saturate due to nonlinear effects until a ~ 3. In other words, 
the r-modes in a strange star may never actually reach the 
truly nonlinear regime. If this is, indeed, the case our spin- 
down model may be better than we initially had any reason 
to expect. 

We now consider the implications of the above results 
for the r-modes as a source of detectable gravitational waves. 
We estimate the g ravita tional- wave strain using [cf. Ander- 
sson & Kokkotas ( ^00 1| )], 

h « 7.54 X 10-^'' aJPliMiARw^^-^ (32) 

where D is the distance to the source (here taken to be in 
the Virgo cluster). Typical results are shown in Figure ^ 
The oscillatory evolution of the r-mode amplitude would, 
obviously, imprint a unique signature on the emerging grav- 
itational waves. In contrast to the neutron star case, where 
the r-mode signal evolves slowly on the cooling timescale, 
the signal from a strange star would be concentrated in rel- 
atively brief "bursts" . Even though one must recognize the 
difficulties in providing a detailed model for this behaviour 
(given the need for a detailed, probably nonlinear, under- 
standing of for example the viscosities) it seems clear that 
observations ought to make it quite easy to distinguish a 
strange star from a neutron star. 

We will not attempt to assess the detectability of the 
early phase of the gravitational-wave signal shown in Fig- 
ure 0. The rapid (model dependent!) variations in the mode 
amplitude implies that this problem requires a more detailed 
signal processing discussion. However, this initial phase last 
only for a few years. Once the star reaches the quasi- 
adiabatic phase where a remains constant over long peri- 
ods in time we can easily estimate the detectability of the 



gravitational waves. To do this we consider the data in Fig- 
ure ^. From the figure which we infer that h ~ 10~^^ and 
/ ~ 900 Hz at t ~ 10** s . Then, knowing that the effective 
amplitude improves roughly as the squareroot of the number 
of detected cycles, i.e. that hes ~ V/^obsft, we find that one 
would need an integration time of fobs ~ 40 years to achieve 
an effective amplitude of order 10~^^ . This is clearly not fea- 
sible, and we can conclude that the late part of the signal 
shown in Figure ^ is unlikely to be detected from sources in 
the Virgo cluster. 

The usual rationale for considering sources in the Virgo 
cluster is that we need a reasonable event rate to make a 
strong case that a certain model leads to detectable gravita- 
tional waves. In the present case, it may be sufficient to focus 
on sources in our own Galaxy. After all, the gravitational- 
wave signal from a strange star lasts considerably longer 
than that from a young neutron star. From Figure H we can 
see that the signal remains above h ~ 10~^^ for more than 
300 years. Taking the corresponding wave frequency to be 
600 Hz (again extrapolated from the figure), we find that the 
effective amplitude would reach 10~^^ after two weeks of in- 
tegration if the source was located at a distance of 10 kpc, 
i.e. at the centre of the Galaxy. Finally combining the ex- 
pected supernova rate on 3/yr/galaxy with a lifetime of the 
gravitational wave signal of 300 years or so, we deduce that 
there could be as many as 10 or so systems active in the 
Galaxy at any given time. At least if all compact stars born 
in supernovae are strange. 

The above results offer exciting possibilities given the 
fact that new large scale interferometers are about to come 
into operation. Our study suggests that these detectors 
should be able to test the idea that "all neutron stars are 
strange stars" in the near future. 

Before we conclude our discussion of young strange 
stars, it is worthwhile making two remarks. Firstly, one 
might wonder whether the r-mode driven spin-down would 
have observational effects in addition to the generated gravi- 
tational waves. It could, for example, leave a distinct imprint 
on the electromagnetic pulsar signal. After all, the r-mode 
phase may last for thousands of years in a young strange 
star born spinning at, or near, the Kepler limit and it is 
possible that the presence of the gravitational-wave torque 
could lead to some anomalous spin-down behaviour. Unfor- 
tunately, one can argue that this is unlikely to be the case, 
unless the star is observed in the first 10 years or so. One 
can easily show that the r-mode spin-down torque is domi- 
nated by the standard contribution from an electromagnetic 
dipole, 

(33) 

whenever 

B > 1.7 X W^^aMiAPIsG (34) 

Taking the spin-period to be 2 ms and assuming a magnetic 
field of 10^^ G (fairly typical for a young pulsar) we find 
that the spin-down will be mainly electromagnetic once the 
mode-amplitude falls below a ~ 10~^. From Figure ^ we can 
see that this will happen after the first ten years or so. In 
other words, we should not necessarily expect the presence 
of a small amplitude r-mode to make an imprint on the 
electromagnetic signal from a young pulsar. 
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Secondly, we note a recent paper by Middleditch et al 



(2000) which provides possible evidence for the existence 



of a 2.14 ms pulsar in the remnant of SN1987A. The ob- 
servations were made over several years between 1992 and 
1996, and the data indicate a surprisingly large, and highly 
variable spin-down rate. It is interesting to note that these 
suggestions agree qualitatively with the results of our sim- 
ple spin-evolution model for a newly born strange star. In 
particular, it is clear that one would expect to find a highly 
variable P during the phase where the r-mode amplitude 
oscillates on the cooling timescale. Should the Middleditch 
observations prove to be reliable, they may thus indicate 
that a strange star was born in SN1987A. 



6 SUMMARY 

In this paper we have studied the r-mode instability in the 
context of strange stars. We have shown that unstable r- 
modes affect strange stars in a way that is quite distinct from 
the neutron star case. For accreting strange stars, the onset 
of r-mode instability does not lead to a thermo-gravitational 
runaway. Instead, the strange star evolves towards a quasi- 
equilibrium on a timescale of about a year. This could ex- 
plain the clustering of spin-frequencies inferred from kHz 
QPO data in Low-mass X-ray binaries. For young strange 
stars we showed that the r-mode driven spin-evolution is 
also distinct from the neutron star case. In a young strange 
star the r-mode undergoes short cycles of instability during 
the first few months. This is followed by a quasi-adiabatic 
evolution where the r-mode remains at a small, roughly con- 
stant, amplitude for perhaps as long as 10^ years. Another 
interesting feature of these evolutions is that the r-modes in 
a strange star never grow to large amplitudes. This could 
prove to be a crucially impo rtant observation since recent 
hydrodvnamical simulations ( Stergioulas fc Font 2000 Lind- 



blom, Tohline fc Vallisneri 200o| ) indicate that the r-modes 
will not be significantly affected by nonlinear effects until at 
much larger amplitudes than those reached in our study. 

The main conclusion of this study is that the r-modes 
in a strange star should emit a persistent gravitational-wave 
signal that ought to be detectable with large-scale interfer- 
ometers given an observation time of weeks to months. If 
detected, these signals would provide unique evidence for 
the existence of strange stars in the Universe, which would 
put useful constraints on the parameters of QCD. 
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APPENDIX: TWO COMPLICATING ISSUES 

In this Appendix we discuss briefly two issues that could 
prove to be relevant for the spin-evolution of a strange star. 



A preliminary study shows that these effects will not alter 
the conclusions drawn in this paper qualitatively, but that 
they must be considered in any attempt to obtain reliable 
quantitative estimates. 



6.1 Deconfinement of the crust 

As mentioned in the Introduction, a strange star is likely 
to have a thin crust of normal nuclear matter. The inner 
edge of this crust is sharply defined by the nuclear drip den- 
sity, since the crust relies on the electromagnetic field to 
suspend it above the strange matter core and save it from 
deconfinement. Neutrons are obviously not affected by the 
magnetic field and will fall into the strange core as soon as 
they begin to drip out of the nuclei. For a given equation of 
state, one can determine the maximum mass of the crust. 



This has been done by Glendenning and Weber (1992). Not 



very surprisingly, the maximum mass dep ends a lso on the 
rotation rate of the star. Yuan and Zhang (1999) have pro- 
vided a simple analytic calculation of the maximum crust 
mass, which agrees reasonably well with the numerical re- 
sults of Glendenning and Weber. For simplicity, we will use 
the analytic result here, i.e. we assume that 



crust, max ■ 



[2.2 + 2.9R^o^M~l{n/flK, 



(35) 



Clearly, this estimate may have impact on a discussion of 
accreting strange stars. 

Strange star accretion can be seen as a two-stage pro- 
cess. First, the accreted matter, of mass say dM, is added 
to the crust. If Mcrust < Mcrust.max the increment dM can 
be added to the crust without any matter migrating into 
the core and thus no significant accretion heating. This is 
particularly relevant since accretion leads to spin-up and an 
increase in Mcrust , max 

Conversely, in the case of a rapidly spinning down star 
(perhaps due to an unstable r-mode), a substantial por- 
tion of the crust will 'fall' into the core on the spindown 
timescale. This may then rapidly heat the core. Consider 
the case of spindown due to a saturated r-mode, for which 



n ■ 



Taking the time-derivative of (pq) we have 



Mc ~ 5.85 X IQ-^'R^o^Mri^ Mq/s 



or, after using (p6[), 

Mc = -1.5 X 10"'' a'^ Rio Mi l PI^ Mq/s 



(36) 



(37) 



(38) 



For large mode amplitudes and a spin-rate near the Kepler 
limit, this would correspond to a huge rate of deconfinement 
of the crust. The effect of this on the temperature of the 
star, and the evolution of the unstable r-mode could clearly 
be significant. 

The above discussion does, however, come with a few 
caveats. The electromagnetically supported gap between the 
crust and core is incredibly thin: Alcock et al. (1986) sug- 
gest 10~^^ m, which is of order 10"'^ atomic diameters, or 
10 nuclear diameters. This means that, unless the crust and 
core execute nearly exactly the same oscillation the two will 
overlap with part of the crust "dissolving" into the core. In 
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the hmit where the crust is infinitely rigid and so does not 
partake in the oscillation, we can imagine the interior oscil- 
lating core dissolving the inner part of the crust, in a shell of 
thickness corresponding to the radial motion of the strange 
core. The crust would then presumably rapidly contract, 
with another shell being dissolved, and so on. Of course, in 
reality the crust is far from rigid. Being much thinner and 
less massive than a neutron star crust, its modulus of rigidity 
is very small. It follows that the crust may participate in a 
fluid oscillation mode to a considerable extent, [as discussed 
in the neutron star case by Levin & Ushomirsky (2000)]. 
We have not attempted to quantify the extent to which the 
crust of a strange star participates in (say) a core r-mode, 
but it is an interesting question that may be worth further 
consideration. 

Finally, the fact that the core-crust boundary corre- 
sponds to the density at which neutrons being to drip out of 
the nuclei implies that possible viscous Ekman layers will be 
much less significant th an in the neutron star case, cf. the 
comments by Madsen (2000). In view of this, we have not 
included effects due to viscous rubbing at the crust interface 
in our analysis. 




logioT(K) 

Figure 5. An illustration of the relevance of the nonlinear contri- 
bution to the strange star bulk viscosity. We show the curves in 
the P — T plane corresponding to various values of the saturation 
amplitude Os, cf. Eqn Also shown is the main instability 

window for a strange star (thick solid line). 



6.2 Nonlinear viscosity saturation 

The possibility that nonlinear contributions to the viscosity 
coefficients may saturate an unstable mode has been sug- 
gested by Reisenegger (private communication). However, 
in the case of a "normal" neutron star fluid the standard 
bulk viscosity is too weak for this effect to be relevant (the 
mode does saturate, but only at unphysical amplitudes). In 
the case of strange stars, the bulk viscosity is significantly 
stronger and thus it could well be that the nonlinear contri- 
butions to the viscosity turn out to be relevant. Incidentally, 
it may be worth pointing out that a similar effect may be 
highly rel evant in the case of strong hyperon induced bulk 



viscosity (lones 2001; Lindblom & Owen 2001) 



We want to investigate the consequences of the fact that 
the low temperature bulk viscosity coefficient can be written 
(cf. eqn. (21) of Madsen ([l992|)) 



, ,2 



1 + 



167r2 \?,^ldVoT 



(39) 



In the estimate used earlier in the paper (eg. to create the 
data displayed in Figure |l|) the Av/vo term was neglected. 
The bulk viscosity calculation could, of course, be gener- 
alised in order to include this compression term but before 
doing this it is useful to estimate whether this contribution 
will ever be relevant. To do this we use 

— - I (40) 

where 8r is the radial component of the r-mode displace- 
ment vector. The two terms in the bracket of ( p9| ) contribute 
equally when 



Sr 



0.07a 



7.5 X 10 ^^4ooniioo29 



(41) 



where we have used the estimate of the height of the r- 
mode at t he sur face of the star provided by Andersson and 
Kokkotas (2001). Parameterising to typical values we can 
use this result to get an estimate of the r-mode amplitude 



at which the nonlinear contribution to the bulk viscosity 
should not be neglected. We then get 



1.07T9//4oomioo (;^) 



(42) 



from which we can deduce that the contribution may be 
important in stars cooler than (say) 10^ K, i.e. in the main 
instability of Figure |l[ 

We can also estimate the amplitude at which the r- 
mode saturates because of the nonlinear contribution to the 
viscosity. To do this (in a suitably simple way that avoids 
volume integration over eigenfunctions etc) we use 



5r 



167r2 \^3^idr R 



: 0.36Q^P_3'miooM4oor| 



and modify the bulk viscosity timescale to 



tt™ = tb\l + 0.36a^P_3 miooMlooTf] 



(43) 



(44) 



Assuming that the mode amplitude is large (such that we 
can neglect the first term in the bracket of the denomina- 
tor), and that the shear viscosity is suitably weak, we arrive 
at a saturation amplitude (by setting a — the evolution 
equation for the mode amplitude) 



19.4^i| 



(45) 



Although this is a rough estimate that may be far from the 
detailed answer it shows that the effect could, in principle, 
be relevant and that nonlinear viscosity saturation should 
be taken seriously for strange stars. 

However, as is clear from Figure ^ a young strange star 
never evolves into the regime where the nonlinear contribu- 
tion to the bulk viscosity dominates. Hence, although it is 
conceptually interesting, we do not expect this effect to play 
a role in saturating an unstable r-mode in a young strange 
star that enters the instability window spinning at or near 
the breakup limit. 



10 N. Andersson, D.I. Jones and K.D. Kokkotas 



REFERENCES 

Alcock C, Farhi E., Olinto A., 1986, Ap. J. 310 261 
Alpar, M.A., 1987, Phys. Rev. Lett. 58 2152 

Andersson, N., Kokkotas, K.D., 2001, Int. J. Mod. Phys. 10 381 
Andersson N., Jones D.I., Kokkotas K.D., Stergioulas N., 2000 

Ap. J. Lett. 534 L75 
Andersson N., Comer G.L., 2001, MNRAS, 328, 1129 
Andersson N., Kokkotas K.D., Schutz B.F., 1999, Ap. J 510 846 
Andersson N., Kokkotas K.D., Stergioulas N., 1999, Ap. J 516 

307 

Bildsten L., 1998 Ap. J. Lett. 501 89 

Bildsten L., Ushomirsky G., 2000, Ap. J. Lett. 529 33 

Colpi M., Miller J.C., 1992, Ap. J. 388 513 

Edwards R.T., van Straten W., Bailes M., 2001, Ap. J. 560 365 
Farhi E., Jaffe R.L., 1979, Phys Rev D 30, 2379 
Glendenning N. K., Weber F., 1992, Ap. J. 400 647 
Ho W.C.G, Lai D., 2000, Ap. J. 543 386 
Jones P.B., 2001, Phys. Rev. Lett. 86 1384 

Lattimer J.M., van Riper K.A., Prakash M., Prakash M., 1994, 

Ap. J. 425 802 
Levin Y., 1999, Ap. J. 517 328 



Levin Y., Ushomirsky, G., 2000, 3.stro-ph/0006028 



Lindblom L., 2001, Neutron Star Pulsations and Instabilities pp. 
257-276 in "Gravitational Waves: A Challenge to Theoretical 
Astrophysics," Ed; V. Ferrari, J.C. M iller, and L. Rczz oUa 
(ICTP, Lecture Notes Series) [preprint [sistro-ph/OlOllSet 

Lindblom L., Owen B.J., Morsink S., 1998, Phys. Rev. Lett. 80, 
4843 

Lindblom L., Owen B.J., Ushomirsky G., 2000, Phys Rev D. 62 
084030 

Lindblom L., Owen B.J., 2001, Effe ct of hvveron bulk viscosity 



on neutron-star r-modes preprint astro-ph/011055S 
Lindblom L., Mendell G., 2000, Phys. Rev. D 61 104003 
Lindblom L., Tohline J.E., Vallisneri M., 2001, Phys. Rev. Lett. 

86 1152 

Madson J., 1992, Phys. Rev. D 46 3290 
Madsen J., 1998a, Phys. Rev. Lett. 81 3311 

Madsen J., 1998b, in "Hadrons in Dense Matter and Hadrosyn- 
thesis". Lecture Notes in Physics, Springer Verlag, Ed. 
J.Cleymans 
Madsen J., 2000, Phys. Rev. Lett. 85 10 
Mendell G., 2001, Phys. Rev. D 64 044009 
Middleditch J., et al. 2000, New Astronomy 5 243 
Owen B.J., Lindblom L., Cutler C, Schutz B.F., Vecchio A., An- 
dersson N., 1998, Phys Rev D, 58, 084020. 
Page D., Prakash M., Lattimer J.M., Steiner A.W., 2000, Phys. 

Rev. Lett. 85 2048 
Papaloizou J., Pringle J.E., 1978, MNRAS 184 501 
Pizzochero P. M., 1991, Phys. Rev. Lett. 66 2425 
RezzoUa L., Lamb F.K., Shapiro S.L., 2000, Ap. J. Lett. 531 139 
Spruit H.C., 1999, Astron. Astrophys. 41 LI 
Stergioulas N., Font J.A., 2001, Phys. Rev. Lett. 86 1148 
Usov, V.V., 2001, Phys. Rev. Lett. 87, 1101 
Wagoner R.V., 1984, Ap. J. 278 345 

Watts A.L., Andersson N. . 2001 The spin evolution of nascent 
neutron stars preprint [istro-ph/0110573| 



Witten E., 1984, Phys. Rev. D 30 272 
Yuan Y. F., J. L. Zhang, 1999, Astron. Astrophys. 344 371 



